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FOREWORD 


This  report,  "A  Digital  Simulation  Program  Describing  the 
Motion  of  an  Aircraft  Undergoing  Engine  Failure  During  its 
Takeoff  Ground  Roll,"  describes  the  primary  actions  of  the 
program  ENGOUT,  and  serves  as  a  User's  Guide  for  the  program. 
This  program  was  written  to  give  stability  and  control  engineers 
a  tool  to  predict  the  ground  minimum  control  speed  of  a  multi¬ 
engined  aircraft.  The  work  resulting  in  this  report  was  begun 
in  June  of  1980  and  completed  in  March  of  1981  by  Mr  Michael  J. 
Miedlar ,  Project  Engineer-,-  Flight  StaBitniy'~arKa"Control  Branch 
(ASD/ENFTC) .  This  report  was  submitted  by  the  author  in  July 
of  1981. 
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SECTION  I 


INTRODUCTION 


The  problem  of  engine  failure  during  takeoff  enters  into 
the  design  phase  of  all  multi-engined  aircraft.  Ground  minimum 
control  speed  (Vmcg)  has  a  direct  impact  on  sizing  the  vertical 
tail,  the  rudder  and  its  deflection  authority,  as  well  as  on 
any  directional  stability  augmentation  system.  The  difficulty 
associated  with  this  problem  arises  from  the  estimation  of 
landing  gear  tire  friction  and  reaction  forces.  These  forces 
depend  on  factors  such  as  runway  surface  condition,  speed,  and 
airplane  weight. 

This  report  is  intended  to  serve  as  an  introduction  to  the 
program  "ENGOUT",  which  solves  this  problem.  This  program, 
written  for  use  on  a  "MODCOMP  CLASSIC"  digital  computer#  generates 
a  time  history  which  describes  the  forces,  moments,  and  motion 
resulting  from  an  engine  out  condition  during  the  takeoff  ground 
roll.  This  program  will  analyze  any  airplane  which  has  a  nose 
wheel  landing  gear  bogie  and  up  to  four  main  gear  bogies. 

"ENGOUT"  will  accept  as  many  as  four  engines,  located  anywhere 
on  the  structure.  Any  combination  of  these  engines  may  fail, 
at  the  same  time,  at  any  speed  during  the  ground  roll.  Crosswind 
may  be  entered  as  a  ramp  or  a  step  function  from  either  side  of 
the  runway.  Center  of  gravity,  gross  weight,  runway  condition, 
and  pilot  reaction  time  are  all  variable. 

The  program,  written  in  the  "MIMIC"  programming  language 
for  the  KC-135  re-engining  project,  has  been  generalized  to 
enable  it  to  analyze  a  wide  variety  of  aircraft  and  takeoff 
conditions.  A  brief  description  of  "MIMIC"  is  given  in  Appendix 
A. 
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SECTION  II 
USING  ENGOUT 


ENGOUT  generates  the  time  history  of  an  aircraft  during  a 
takeoff  ground  roll.  The  user  may  choose  to  fail  an  engine  at 
any  speed  during  the  ground  roll,  and  the  program  will  predict 
the  trajectory  of  the  aircraft.  The  plane  yaws  in  the  direction 
of  the  failed  engine,  and  a  simulated  pilot  enters  control  to 
return  the  aircraft  to  its  original  flight  path.  ENGOUT  follows 
the  aircraft  to  rotation  speed  if  no  failure  occurs,  or  to  the 
maximum  point  of  deviation  following  an  engine  failure.  Although 
able  to  analyze  a  routine  takeoff  roll,  ENGOUT  is  primarily 
intended  for  use  as  a  tool  in  determining  Vmcq»  the  ground  minimum 
control  speed.  y 

The  ground  roll  does  not  need  to  start  at  zero  velocity, 

ENGOUT  has  the  capability  to  begin  at  any  speed  less  than  or 
equal  to  the  speed  at  which  failure  occurs.  Any  speed  of  cross- 
wind  may  be  entered  as  a  ramp  or  step  function  at  any  value  of 
forward  speed,  from  either  direction.  Thrust  is  a  function  of 
speed,  and  the  program  provides  the  capability  to  simulate  non- 
linearities  in  this  function.  The  program  also  has  capability 
to  describe  the  dependency  of  maximum  rudder  deflection  on 
aircraft  speed.  The  program  has  enough  flexibility  to  account 
for  a  great  variety  of  phenomena. 

ENGOUT  can  analyze  almost  any  conceivable  takeoff  condition. 
Variable  parameters  include  crosswind,  center  of  gravity,  gross 
weight,  runway  surface  condition,  failed  engine  spin  down  time 
constant,  or  any  other  aircraft  parameter  of  interest.  Written 
in  the  hybrid  language  MIMIC,  this  program  manipulates  forces 
and  moments  to  calculate  linear  motions  in  the  X  and  Y  inertial 
runway  directions  and  angular  motions  in  the  roll  and  yaw  planes. 
Using  MIMIC's  implicit  integration  routine,  it  integrates  these 
accelerations  to  determine  the  velocities  in  the  X  and  Y  direc¬ 
tions,  the  roll  rate,  and  the  yaw  rate.  Integrating  again,  ENGOUT 
calculates  the  X  and  Y  distances,  the  roll  angle  4>,  and  the  yaw 
angle  ij>.  It  takes  into  account  all  aerodynamic  forces  and 
moments,  as  well  as  all  reactions  and  friction  forces  acting  on 
the  landing  gears.  While  on  the  ground,  an  aircraft  does  not 
translate  in  the  Z  direction,  nor  does  it  pitch.  ENGOUT  solves 
the  four  degree  of  freedom  problem  of  an  aircraft  during  its 
ground  roll. 

ENGOUT  provides  input  and  output  variable  dictionaries.  The 
user  enters  data  for  the  problem  of  interest.  This  is  done  by 
the  techniques  used  in  MIMIC  as  described  in  appendix  A.  The 
data  must  be  entered  in  the  MIMIC  format,  and  in  a  particular 
set  of  units.  All  lengths  must  be  in  feet.  Enter  all  areas  in 
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ft2,  all  angles  in  degrees,  and  all  inertias  in  slug/ft2.  Mass 
should  be  in  slugs,  density  in  slugs/ft3,  time  in  seconds,  and 
velocity  in  knots.  ENGOUT  works  with  all  stability  and  control 
derivatives  in  1/rad.  The  program  works  with  data  as  though 
entered  in  these  units,  so  the  user  will  need  to  ensure  the 
accuracy  of  the  dimensions  of  his  data. 

When  obtaining  geometric  data  for  the  aircraft  of  interest, 
use  a  three  view  drawing  which  describes  the  vehicle  in  fuselage 
stations,  butt  lines,  and  water  lines.  The  program  manipulates 
these  terms  to  obtain  the  necessary  moment  arms.  One  term  will 
present  some  difficulty,  the  fuselage  station  of  each  engine. 

This  term  describes  the  point  at  which  the  vertical  component  of 
an  engine's  thrust  acts  in  the  pitch  plane.  The  engines  attach 
to  the  structure  in  some  way,  and  the  vertical  thrust  component 
transmits  a  pitching  moment  thru  this  structure.  This  moment 
arm  defines  the  fuselage  station  of  the  engine.  All  fuselage 
stations,  water  lines,  and  butt  lines  should  be  entered  in  feet. 

The  program  looks  for  data  on  four  main  gear  bogies  and  four 
engines.  If  the  airplane  being  analyzed  does  not  have  these 
components,  enter  a  zero  for  any  parameter  to  do  with  the 
component.  The  program  has  logical  controls  to  account  for  this 
option. 

ENGOUT  returns  information  in  a  particular  set  of  units  also. 
All  angular  measures  output  in  degrees,  and  all  linear  measures 
in  feet.  All  forces  will  have  units  of  pounds.  Logical  terms 
are  non-dimensional.  Any  output  variables  pertaining  to  compo¬ 
nents  not  present  on  the  aircraft  will  have  values  of  zero.  The 
only  exception  to  this  is  logical  control  variables  for  non-present 
components.  Some  of  these  will  have  a  true  value  of  one,  and  some 
the  false  value  zero.  The  program  itself  determines  the  value, 
which  will  depend  upon  the  use  of  the  variable.  The  variables 
listed  in  the  output  dictionary  are  recommended  to  the  user  as 
completely  able  to  describe  the  trajectory  of  the  aircraft.  Of 
course  the  actual  output  terms  are  at  the  discretion  of  the  user. 

ENGOUT  works  like  a  black  box.  The  user  puts  in  data  and 
gets  back  information  describing  the  trajectory  of  the  aircraft. 

The  following  sections  of  this  report  describe  the  main  actions 
of  the  program.  It  looks  inside  of  the  black  box.  Appendix  B 
lists  the  program,  and  shows  a  complete  run  for  a  certain  set  of 
data. 
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SECTION  III 


GROUND  MINIMUM  CONTROL  SPEED 


From  being  initially  at  rest,  an  aircraft  will  accelerate 
down  the  runway  under  takeoff  power.  After  reaching  Vs,  the 
stall  speed,  the  aircraft  is  capable  of  flight,  but  the  pilot 
will  not  begin  rotation  until  the  aircraft  has  accelerated  past 
the  critical  engine  failure  speed,  Vj_.  The  critical  engine 
failure  speed  is  the  maximum  speed  at  which  the  pilot  can  control 
the  aircraft  after  the  failure  of  a  critical  engine  and  either 
takeoff  or  bring  the  aircraft  to  a  stop  within  the  runway  availa¬ 
ble.  Above  this  speed  the  pilot  will  rotate  the  aircraft  and  lift 
off. 


Safe  takeoff  procedure  requires  that  the  pilot  knows  the 
critical  engine  failure  speed.  However,  a  more  crucial  parameter 
to  the  design  engineer  is  the  ground  minimum  control  speed,  Vmcg. 
This  is  the  minimum  speed  at  which,  if  a  critical  engine  fails, 
the  pilot  should  be  able  to  regain  control  of  the  aircraft  within 
a  certain  specified  deviation  from  the  runway  centerline,  and 
continue  the  takeoff.  Normal  design  procedures  require  Vmc  to 
be  less  than  or  equal  to  Vi. 

Naturally,  a  low  vmcg  is  very  desirable,  so  a  manufacturer 
will  do  all  he  can  to  minimize  this  speed.  By  choosing  a  favorable 
takeoff  condition,  the  minimum  control  speed  may  indeed  appear  very 
low.  But  a  seemingly  inconsequential  change  in  one  of  the  factors 
involved  in  the  analysis  may  have  significant  impact  on  its  results. 

An  engine  failure  aggravates  the  already  complicated  problem 
of  takeoff  analysis.  The  reaction  forces  on  the  landing  gear  tires 
must  be  calculated  accurately  to  ensure  a  proper  result.  Entities 
such  as  airplane  gross  weight,  center  of  gravity  position,  runway 
surface  condition,  and  even  crosswind  velocity  all  influence  these 
reactions.  Each  of  these  factors  must  enter  into  any  analysis  of 
the  phenomenon. 

By  far,  the  time  for  pilot  reaction  has  more  impact  than  any 
other  factor.  The  faster  he  reacts  to  the  failure,  the  more 
controlled  the  motion  will  seem,  and  in  fact  will  be.  Controversy 
rages  around  this  time  delay,  and  no  standard  has  yet  been  univer¬ 
sally  agreed  upon.  A  small  change  in  reaction  time  may  result  in 
a  large  change  in  aircraft  motion,  and  give  the  appearance  of 
better  performance. 

One  manufacturer  tries  to  get  around  this  problem  by  installing 
a  special  sort  of  directional  stability  augmentation  system.  In 
effect,  the  aircraft  itself  reacts  to  the  failure,  and  puts  in 
some  control.  The  re-enginer  KC-135  aircraft  is  now  equipped  with 
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engine  failure  command  augmentation,  and  this  type  control  input 
must  also  be  accounted  for  in  an  analysis  of  Vmcg,  as  well  as  the 
action  of  a  standard  directional  augmentation  system  as  installed 
on  most  large  aircraft. 

These  factors  all  interconnect  to  make  the  problem  a  most 
difficult  one  to  solve.  The  input  parameter  of  primary  interest 
is  the  speed  at  which  the  engine  fails,  and  the  output  parameter 
of  especial  interest  is  the  deviation  from  the  runway  centerline, 
since  these  are  what  define  the  minimum  control  speed.  VmCg  is 
the  speed  at  which  a  critical  engine  failure  can  occur,  and  after 
a  reasonable  time  delay,  the  pilot  can  control  the  aircraft  before 
it  has  strayed  a  specified  distance  from  the  runway  centerline. 

A  critical  engine  is  defined  as  that  engine  whose  failure  would 
cause  the  largest  yawing  moment,  that  is,  the  most  outboard  engine, 
since  engines  are  mounted  symmetrically  in  the  yaw  plane. 

The  calculation  of  vmca  involves  these  factors,  and  many 
others.  To  solve  this  problem,  the  forces  and  moments  acting  on 
the  aircraft  have  to  be  computed,  and  integrated  for  its  motion. 

The  ground  physically  restricts  the  vertical  motion  of  the  landing 
gear  wheels,  and  the  gear  bogies  are  assumed  to  be  rigid.  The 
pitch  attitude  is  held  constant.  Therefore,  the  aircraft  has  four 
degrees  of  freedom,  and  these  motions  all  contribute  to  the  total 
motion  of  the  stricken  vehicle. 
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SECTION  IV 


CALCULATION  OF  TIRE  NORMAL  WHEEL  FORCES  (Fz) 


During  any  ground  roll,  the  forces  acting  on  the  landing 
gear  tires  of  an  aircraft  greatly  affect  its  trajectory  down 
the  runway.  Normal  forces  support  its  weight,  rolling  friction 
acts  as  a  drag,  and  side  forces  oppose  translation  to  either 
side.  These  frictions  and  reactions  also  oppose  any  rotation 
of  the  aircraft,  which  directly  influences  the  aircraft  motion. 

Of  primary  importance  are  the  gear  reaction  forces,  the 
normal  forces  which  support  the  weight.  These  reactions  enter 
into  the  calculation  of  the  friction  forces,  and  prevent  an 
aircraft  from  pitching  during  a  ground  roll.  These  reactions 
depend  on  center  of  gravity  location,  since  all  forces  and 
moments  can  be  visualized  as  acting  at  the  center  of  gravity. 
Proper  manipulation  of  these  forces  and  moments  leads  to  the 
calculation  of  the  gear  reaction. 

ENGOUT  uses  a  basic  equation  commonly  known  as  the  three 
moment  equation,  to  calculate  these  reactions.  It  describes 
the  relationship  between  the  internal  bending  moments  at  three 
consecutive  points  of  support  on  a  beam  continuous  over  three 
or  more  supports,  as  shown  in  figure  1,  under  any  external 
loading  condition.  This  equation  is  used  in  conjunction  with 
statics  to  determine  reaction  forces  for  a  statically  indeter¬ 
minate  beam. 

Assuming  that  the  landing  gears  act  in  the  pitch  plane  as 
a  series  of  simple  support,  the  use  of  this  equation  is  valid. 

In  the  pitch  plane,  gears  with  the  same  fuselage  station  have 
the  same  moment  arms,  and  act  as  one.  In  general,  this  will 
reduce  the  number  of  supports  to  two  if  the  aircraft  has  one 
pair  of  main  gear  bogies,  or  three  if  it  has  two  pairs  of  main 
gear  bogies.  Because  gear  bogies  with  the  same  fuselage  station 
may  be  combined  into  one  support,  aircraft  are  built  with  landing 
gear  bogies  symmetrically  placed  about  the  fuselage  centerline. 

Applying  this  assumption,  an  aircraft  resembles  a  simply 
supported  rigid  beam,  making  it  possible  to  calculate  the  normal 
reactions.  Two  equations  available  from  statics  are  the  summation 
of  moments  in  the  pitch  plane  and  the  summation  of  forces  in  the 
vertical  direction.  If  the  plane  has  two  points  of  support, 
specifically  a  nose  gear  logic  and  two  main  gear  bogies  acting 
together,  these  are  sufficient  to  calculate  the  reactions. 

However,  three  points  of  support,  a  nose  gear  bogie  and  two  sets 
of  main  gear  bogies,  require  an  additional  equation,  the  three 
moment  equation. 
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Consider  the  beam  continuous  over  three  supports  shown  in 
figure  2.  Ma,  Mg,  and  Me  represent  the  bending  moment  at  each 
of  the  supports  due  to  any  potential  loading  of  the  entire  beam. 
Fi  and  Fj  represent  any  potential  loading  condition  acting  on 
the  partial  spans  and  L2.  This  method  treats  each  length  of 
beam  Li  and  L2  independently  from  the  rest  of  the  beam.  Each 
length  of  beam  has  an  associated  bending  moment  diagram  generated 
by  the  particular  loading  condition  acting  on  it  alone.  Suppose 
that  the  bending  moment  diagrams  have  an  abritrary  shape  as  shown 
in  figure  3. 

Ai  and  Aj  represent  the  areas  associated  with  the  bending 
moment  diagrams  of  each  segment,  a^  and  Ej  are  lengths  from 
the  indicated  ends  to  the  geometric  centroid  of  each  area. 

Based  on  these  parameters,  the  three  moment  equation  as  listed 
in  reference  5  is  given  by  equation  (1) . 

Ma  +•  2  Mjj  (L^  +  La)  +  Mg  L2  = 

m  n  (1) 

—  Z  6  Aj_  ai  —  Z  6  A-j  b-j 

1=1  ~tr~  j=1 


"m"  and  "n"  are  the  number  of  distinct  loads  acting  on  each 
length  of  beam.  For  a  beam  with  more  than  three  points  of 
support,  this  equation  can  be  repeated  as  many  times  as  necessary 
to  include  the  bending  moments  a  teach  support  along  the  entire 
beam,  taking  three  adjacent  supports  each  time.  These  equations 
can  be  solved  simultaneously  to  determine  the  bending  moments  at 
each  support.  These  moments,  along  with  the  equations  for  summa¬ 
tion  of  forces  and  moments,  can  be  used  to  calculate  the  normal 
reaction  forces  at  each  support. 

Depending  on  the  loading  condition,  the  three  moment  equation 
takes  on  special  forms .  Suppose  that  the  loads  consist  of  concen¬ 
trated  vertical  forces  as  shown  in  figure  4 .  Reference  5  shows 
the  derivation  for  this  special  form  of  the  three  moment  equation 
for  this  loading  condition  as  given  by  equation  (2) . 


Ma  +  2  Mjj  (L^  +  L2)  +  Mg  L2  ~ 

®  2  2  H  ,  2  o 

-  r  pi  aj  (Lx  -  aA)  -  I  pa  bj  (L2  -  bj)  (2) 

i=l  j*1  L2 
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oncentrated  Forces 


The  summation  sign  includes  all  concentrated  loads  on  each 
length  of  beam,  a^  and  bj  are  not  measured  to  the  centroid 
of  the  bending  moment  diagrams  but  from  the  indicated  end  to 
the  point  of  action  of  the  load. 

Suppose  that  a  beam  is  loaded  by  distributed  vertical 
forces  as  shown  in  figure  5.  and  W2  are  the  intensities 
acting  on  each  span.  Reference  5  derives  this  special  form 
of  the  three  moment  equation: 

Ma  Lx  +  2  Mb  (LX  +  L2)  +  Mc  L2  - 


m 

Z 

i=l 


Wi  L1 


n  W4  1*0 

Z  _J _ £ 

j*l  4 


(3) 


The  summation  includes  all  distributed  forces  acting  on  each 
span . 


Figure  6  shows  a  beam  loaded  by  pure  moments  and  Mj . 

The  special  form  of  the  three  moment  equation  for  this  loading 
was  derived  by  the  author  in  the  same  way  that  reference  5  did 
for  the  previous  forms,  and  is  given  by  equation  (4|)  . 

Ma  Lj_  +  2  Mb  (L^  +  I*2)  4-  I»2  * 


m 


-  *  Hi  (Li  -  3  ai  Lx  +  4  a*) 

i=l  T  2 
L1 


(4) 


l  ^  (L2  -  3  bi  L2  +  4  bj) 
j=1  L2 


Once  again,  the  summation  includes  all  pure  moments  acting  on 
each  span.  An  aerodynamic  pitching  moment  is  an  example  of  a 
pure  moment. 

These  three  special  forms  can  be  combined  to  include  as 
many  loading  conditions  as  necessary.  This  combined  special 
form  of  the  three  moment  equation,  taking  into  account  all 
concentrated  vertical  forces,  distributed  vertical  forces,  and 
pure  pitching  moments  is  given  in  equation  (5) . 
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+  2Mh  (Lx  +  I^)  +  McL2  =  (-SPiai  (L^  -  a*)-  Z  Pj  bj  (I^  -  bj) 

L  L2 

+  (-  Z  Wi  L1  -  Z  Wi  h2) 

^  4  4  '* 

3  2  3  (5) 

+  -  Z  **i  (Lx  -  3  ai  1^  +  4  a^ 

2 

Ll 

-  E  *Jl  (L^  ~  3  b[  Iq  +  4  bi) 

2 

L2 

ENGOUT  manipulates  t ha  expanded  form  of  the  three  moment 
equation  to  compute  the  normal  reaction  forces  acting  on  an 
airplane's  landing  gear.  If  the  aircraft  has  only  two  points 
of  support,  each  term  reduces  to  zero,  and  the  equation  loses 
all  significance;  the  program  uses  a  simple  summation  of  forces 
and  moments  to  calculate  reaction  forces.  However,  for  three 
points  of  support,  ENGOUT  calculates  the  appropriate  terms  for 
each  given  loading  condition.  Unnecessary  terms  are  reduced  to 
zero. 

Statics  simplifies  this  process  since  all  forces  can  be 
replaced  by  an  equivalent  force-moment  system  acting  at  a 
different  point.  Lift,  weight,  and  drag,  actually  distributed 
forces,  can  be  transformed  into  three  single  concentrated  forces 
and  one  combined  pitching  moment  acting  at  the  center  of  gravity. 
Thrust,  split  into  its  vertical  and  horizontal  components,  converts 
to  a  set  of  concentrated  vertical  loads  and  pure  pitching  moments. 
The  inertial  reaction  force  also  acts  as  a  pure  pitching  moment  at 
the  center  of  gravity. 

ENGOUT  manipulates  these  forces  and  moments  into  convenient 
forms.  Ma,  the  bending  moment  at  the  nose  gear,  is  assumed  to 
be  zero.  This  is  valid  since,  having  concentrated  the  weight  at 
the  center  of  gravity,  no  forces  or  moments  act  forward  of  the 
nose  gear. 

"ENGOUT"  calculates  M^,  the  bending  moments  at  the  aft  main 
gear,  by  equating  it  to  the  summation  of  moments  aft  of  that  set 
of  gears.  It  then  combines  M^.  with  the  load  terms  from  the  right 
side  of  equation  (5)  to  determine  M3,  the  bending  moment  at  the 
forward  main  gear. 

Having  calculated  the  bending  moments  at  the  three  gear 
positions,  it  is  a  simple  matter  to  calculate  the  reaction 
forces.  Summation  of  bending  moments  forward  of  the  main 
forward  gear  furnishes  the  nose  gear  force.  Summation  of 
bending  moments  forward  of  the  aft  main  gear  produces  the 
reaction  at  the  main  forward  gear  position,  and  summing  the 
vertical  forces  gives  the  main  aft  gear  reaction  force. 
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Although  forces  at  the  same  fuselage  station  act  together 
in  pitch,  they  will  not  do  so  in  roll  since  they  will  be  located 
at  different  butt  lines.  The  roll  moment  arms  will  be  different. 
This  requires  the  division  of  the  reaction  force  between  the 
separate  elements  of  the  support. 

An  airplane's  main  gear  supports  consist  of  two  bogies, 
symmetrically  positioned  below  the  fuselage  centerline.  The 
nose  gear  will  have  only  one  bogie.  This  splitting  of  the 
reaction  force  is  accomplished  by  estimating  equivalent  spring 
force  and  damping  coefficients  for  each  pair  of  bogies.  These 
components  will  act  to  oppose  a  rolling  motion  while  the  wheels 
are  in  contact  with  the  ground.  Also,  the  center  of  gravity  may 
not  necessarily  be  positioned  on  the  fuselage  centerline. 

Letting  "FZMA"  represent  the  total  force  on  the  main  aft  support 
position,  the  reactions  on  each  bogie  is  given  by  equation  (6) . 


F7 

ZMAR 


BLCG  ~  BLMAL  *  FZma  + 
BLMAR  -  CLMFR 


K  *  (BLMAR  -  BLCG)  *  PHI  +  B  (BLMAR  -  BLCG)  *  P 


FZMAL  FzMA  "  FzMAR 


(6) 


"BLCG"  represents  the  butt  line  of  the  center  of  gravity 
position,  "BLMAL"  and  "BLMAR"  the  butt  lines  of  the  main  aft 
right  and  left  gears,  "K"  the  equivalent  spring  force  coeffi¬ 
cient,  "B"  the  equivalent  deunping  coefficient,  "PHI"  the  roll 
angle,  and  "P"  the  roll  rate.  The  corresponding  terms  yield 
the  separate  reactions  for  each  main  forward  gear  bogie. 

ENGOUT  will  perform  these  actions,  thus  delivering  an 
entire  set  of  landing  gear  reactions.  The  equations  necessary 
for  this  procedure  have  been  sufficiently  generalized  to  enable 
the  user  to  analyze  any  aircraft  with  a  nose  gear  and  one  or  two 
sets  of  main  gear  bogies. 
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SECTION  V 


BOLLING  AND  SIDE  FRICTION  FORCES 


During  ground  roll,  an  aircraft's  low  speed  greatly  reduces 
its  aerodynamic  effectiveness,  and  magnifies  the  effect  of 
friction  forces  acting  on  the  landing  gear  tires.  Accurate 
prediction  of  its  trajectory  requires  the  proper  calculation  of 
these  forces.  The  nature  of  friction  makes  an  exact  analysis 
impossible  without  detailed  knowledge  of  the  properties  of  the 
tire  in  question.  However,  it  is  possible  to  give  a  good  approxi¬ 
mation  to  these  forces  by  making  certain  assumptions. 

Rolling  friction  (Fy)  presents  little  difficulty.  Fully 
inflated  tires  deform  under  loading  in  a  characteristic  manner, 
and  a  rolling  friction  coefficient  of  0.025  will  usually  be 
sufficient  to  describe  the  phenomenon.  The  rolling  friction 
force  is  the  product  of  the  normal  force  and  the  rolling  friction 
coefficient. 


FX  =  VR  *  Fz  (7) 

The  rolling  friction  acts  at  the  bottom  of  the  tire,  in  the 
direction  opposite  to  that  of  the  roll. 

Side  friction  force  (Fy)  does  not  lend  itself  to  solution 
quite  so  easily.  This  force  depends  heavily  on  runway  surface 
condition,  tire  yaw  angle,  ground  speed,  and  the  normal  force. 
The  tires  also  react  against  the  side  forces  acting  on  the 
aircraft,  and  if  the  tire  skids,  these  equations  are  invalid. 
The  problems  all  depend  on  the  tire  being  used. 

However,  some  of  these  problems  can  be  disregarded.  The 
reliance  on  tire  yaw  angle  can  be  approximated  as  linear  if 
the  yaw  angle  remains  sufficiently  small,  below  5  degrees. 

Only  test  data  can  provide  totally  accurate  answers  in  this 
region,  but  a  linear  approximation  is  a  good  one. 

Dependency  on  tire  normal  force  can  be  ignored  if  the  yaw 
angle  is  small  enough.  This  is  the  cause  of  the  approximate 
linearity  at  low  yaw  angles. 

As  long  as  a  tire  does  not  skid,  it  will  react  against 
the  side  forces  acting  on  the  aircraft.  It  will  skid  if  the 
side  forces  exceed  the  maximum  static  friction  force.  This 
maximum  is  given  by  the  product  of  the  static  coefficient  of 
friction  and  the  normal  force  acting  on  the  tire.  If  the  tire 
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skids,  the  reaction  force  is  zero  and  the  kinetic  friction 
force  takes  over.  This  force  is  obtained  by  multiplying  the 
kinetic  friction  coefficient  with  the  normal  force.  Reference  6 
presents  a  mathematical  approach  to  landing  gear  ground  reactions 
for  a  general  aircraft,  and  includes  these  equations  to  calculate 
static  and  kinetic  friction  coefficients. 

PK  =  -  0.0053  *  VG  (8) 

pc  -  0.0392  *  RCR  -  0.102  (9) 

“MAX 

These  equations  take  into  account  the  dependency  on  the 
runway  condition  reading  (RCR) ,  which  measures  the  runway 
condition  from  four  for  an  ice  covered  runway  to  23  for  a  dry 
hard  one,  and  VG,  the  ground  speed  in  knots.  Each  coefficient 
has  a  lower  limit  of  0.044. 

ENGOUT  calculates  the  side  forces  acting  on  a  tire  using 
these  assumptions  and  equations.  It  has  logical  control  to 
determine  when  the  side  forces  acting  on  the  airplane  become 
great  enough  to  skid  the  tires.  In  this  case,  the  side  friction 
force  is  set  equal  to  the  kinetic  friction  force. 

The  friction  acting  on  an  unskidded  tire  consists  of  two 
parts,  a  reaction  due  to  the  side  forces  on  the  aircraft  and 
scrubbing  forces  arising  from  the  yaw  angle  of  the  tire.  Each 
tire  reacts  the  same  fraction  of  total  aircraft  normal  force 
and  total  aircraft  side  force.  The  total  side  force  includes 
aerodynamic  side  force  and  centrifugal  force  rising  from  yaw 
velocity. 

The  scrubbing  friction  force  is  the  product  of  the  tire 
normal  force  and  the  scrubbing  friction  coefficient.  The 
scrubbing  coefficient  is  calculated  by  multiplying  the  kinetic 
friction  coefficient  by  a  factor  which  takes  into  account  the 
approximate  linear  dependency  of  friction  with  tire  yaw  angle. 
This  factor  is  determined  by  dividing  the  tire  yaw  angle  by 
eight,  the  maximum  angle  for  which  the  linear  profile  is  valid. 
The  tire  yaw  angle  is  given  by: 

♦tire  “  TAN~1  §  +  LX  *  R/u  (10) 

"V"  represents  linear  velocity  in  the  Y-directionj  "U"  linear 
velocity  in  the  X-direction;  "Lx"  the  distance  in  the  X-direction 
from  the  center  of  gravity  to  the  center  of  the  tire,  positive  if 
the  tire  is  forward  of  the  center  of  gravity,  negative  if  aft; 
and  "R"  the  yaw  velocity  in  radians  per  second. 
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These  assumptions  and  equations  give  a  good  approximation 
of  the  tire  friction  forces.  As  the  aircraft  approaches  liftoff 
speed,  the  effect  of  ground  friction  diminishes  greatly.  However, 
below  ground  minimum  control  speed,  these  forces  have  an  important 
impact  on  an  aircraft's  motion. 
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SECTION  VI 
DIRECTIONAL  CONTROL 


During  a  takeoff  ground  roll,  a  pilot  holds  his  wings 
level  and  keeps  the  aircraft  on  the  runway  center  line  by 
making  small  corrections  to  his  lateral-directional  control 
positions.  These  controls  consist  of  ailerons  and/or  spoilers, 
considered  together  in  wheel/stick  position,  rudder,  and  nose 
wheel  steering  controlled  either  thru  the  rudder  pedals  or  with 
a  tiller.  Basically,  the  pilot  continually  retrims  the  aircraft 
directionally  while  on  the  ground. 

When  an  engine  fails,  its  thrust  spins  down  to  zero  expo¬ 
nentially  and  wind  milling  of  the  propeller  or  fan  adds  drag 
which  rises  exponentially  to  a  steady  state  value .  These  expo¬ 
nentials  have  the  same  time  constant.  This  results  in  a  thrust 
imbalance  in  the  yaw  plane,  causing  the  aircraft  to  nose  around 
in  the  direction  of  the  failure . 

An  aircraft  equipped  with  a  stability  augmentation  system 
senses  instantly  the  motion  caused  by  the  failure  and  puts  in 
rudder  deflection  to  oppose  the  yaw  motion.  The  amount  put  in 
and  how  fast  it  goes  in  depend  on  the  system  installed,  but  this 
sort  of  system  will  put  in  some  control  with  no  recognition  or 
reaction  time. 

A  pilot  needs  some  finite  time  to  recognize,  and  react  to, 
the  situation.  He  then  commands  the  full  directional  control 
available.  He  continues  to  trim  the  aircraft  laterally,  but 
puts  in  full  rudder  and  full  nose  wheel  steering.  He  holds  full 
control  in  until  he  once  again  aligns  the  aircraft  along  the 
runway . 

The  authority  and  effectiveness  of  these  control  systems 
all  depend  on  the  aircraft  being  analyzed,  but  ENGOUT  expands 
and  generalizes  some  factor  to  cover  all  aircraft.  Totally 
aircraft  dependent  features,  such  as  an  augmentation  system, 
can  be  inserted  later  by  the  user. 

The  pre-failure,  or  trimming,  part  of  this  control  is 
simple.  Laterally,  the  pilot  continually  enters  sufficient 
wheel  deflection  to  keep  the  wings  level .  The  block  diagram 
of  figure  1  describes  this  action. 

The  gains  Kp  and  K.  depend  on  the  aircraft  aerodynamics, 
but  values  of  15  for  each  should  be  close.  ENGOUT  models  this 
block  diagram  to  laterally  simulate  a  pilot.  This  "pilot" 
senses  and  responds  to  the  rolling  motion  experienced  during 
the  ground  roll.  P  and  $  influence  P  by  affecting  the  tire 
forces.  Figure  7  does  not  reflect  this,  since  it  only  describes 
the  pilot  lateral  control  action,  valid  for  the  entire  ground  roll. 
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Before  failure,  the  directional  trim  "pilot"  works  the  same 
way,  but  instead  of  "reacting  to"  the  angular  rate  and  angle 
and  entering  the  appropriate  rudder  and  nose  steering  control, 
ENGOUT's  "pilot"  reacts  to  the  side  velocity  and  displacement. 

This  block  diagram  is  shown  in  figure  8 . 

As  before,  the  aircraft  configuration  determines  the  values 
of  Ky  and  Ky.  represents  a  transfer  function  for  converting 

directional  control  to  a  side  acceleration.  ENGOUT  acts  in  such 
a  way  as  to  provide  this  transfer  function.  This  block  diagram 
describes  the  "pilot"  rudder  and  nose  wheel  input  prior  to  an 
engine  failure.  If  the  "pilot"  deflects  both  rudder  and  nose 
wheel  by  moving  the  rudder  pedals,  6ns  will  equal  6r  times  a 
factor  equal  to  the  maximum  nose  steering  angle  available  divided 
by  the  maximum  rudder  deflection  available.  If  a  tiller  controls 
the  nose  wheel  angle  deflection,  each  system  will  need  its  own 
value  of  Ky  and  Ky. 

After  losing  an  engine,  the  pilot  still  holds  his  wings  level, 
and  the  lateral  control  block  diagram  requires  no  additions  or 
changes.  However,  after  recognizing  the  situation,  the  pilot  puts 
in  full  directional  control  in  opposition  to  the  yaw  motion  of  the 
aircraft  created  by  thrust  imbalance.  Once  again,  the  same  block 
diagram  will  describe  both  rudder  and  nose  wheel  steering  systems, 
as  shown  in  figure  9 .  If  both  rudder  and  nose  wheel  are  controlled 
thru  the  rudder  pedals,  the  steering  angle  may  equal  the  rudder 
angle  multiplied  by  the  ratio  of  maximum  nose  wheel  input  to 
maximum  rudder  delfection.  This  will  result  in  maximum  nose 
steering  occur ing  at  the  same  time  as  maximum  rudder  input.  If 
desired,  the  program  will  put  in  maximum  nose  steering  at  a 
different  time  than  maximum  rudder. 

An  augmented  aircraft  will  have  a  lower  value  of  minimum 
control  speed  than  one  without  such  a  system.  The  augmentation 
system  puts  in  rudder  to  oppose  a  yawing  motion  instantly, 
decreasing  the  effect  of  the  delay  between  engine  failure  and 
pilot  reaction.  It  is  usually  designed  with  some  sort  of  washout 
to  ensure  that  it  does  not  oppose  pilot  actions.  The  nature  of 
those  systems,  generally  referred  to  as  yaw  dampers,  is  not  at 
all  general  in  nature,  and  will  be  different  for  each  aircraft. 
ENGOUT  contains  a  simple  system  which  equates  yaw  damper  input 
to  a  gain  times  yaw  rate  plus  a  gain  times  yaw  angle  plus  a  gain 
times  sideslip  angle.  If  a  different  or  a  more  complex  system  is 
needed,  the  necessary  terms  will  have  to  be  inserted  by  the  user. 

An  accurate  engine  failure  simulation  cannot  be  accomplished 
without  proper  simulation  of  the  control  system  of  the  aircraft. 
Unlike  the  similarity  of  the  landing  gears  for  most  aircraft,  the 
control  system  for  each  type  of  vehicle  is  unique.  Especial  care 
should  be  taken  in  modeling  these  systems. 
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APPENDIX  A  -  MIMIC 


The  program  ENGOUT  uses  the  hybrid  language  MIMIC,  a 
language  which  combines  the  speed  of  a  digital  computer  with 
the  programming  ease  of  an  analog  system.  It  eliminates  the 
problems  of  serial  programming  required  by  a  digital  computer 
and  also  the  time  and  amplitude  scaling  usually  needed  for 
analog  operation.  It  uses  a  centralized  fourth  order  variable 
step  Runga-Kutta  routine  for  all  integrations,  making  the  solving 
of  differential  equations  easy.  It  was  developed  to  provide  a 
simple  method  to  solve  equations  on  a  digital  computer.  In 
essence,  the  language  forces  a  digital  computer  to  act  as,  or 
mimic,  an  analog  computer.  MIMIC  can  best  be  described  as  a 
digital-analog  (or  continuous  system)  simulator. 

One  of  MIMIC'S  more  attractive  features  is  its  parallelism, 
which  it  gets  from  its  analog  nature.  The  programming  statements 
need  not  be  in  any  particular  order.  In  this,  it  differs  from 
such  conventional  digital  languages  as  FORTRAN  which  relay  on  the 
order  in  which  statements  execute.  This  allows  MIMIC  programmers 
to  more  closely  match  the  simulations  to  the  real  world  problems 
of  interest. 

MIMIC  uses  a  fixed  format  for  its  statements.  The  result 
always  begins  in  column  10,  and  expressions  begin  in  column  19. 
Columns  2-7  are  reserved  for  logical  variables  which  control 
which  statements  to  execute.  The  FORTRAN  terms  integer  and  real 
number  have  no  meaning  in  MIMIC;  all  numbers  contain  decimal  points. 
Variable  names  can  be  one  to  six  alpha  numeric  characters,  and 
unlike  FORTRAN  may  begin  with  a  number.  The  term  "10."  represents 
the  number  ten,  while  "10"  would  be  interpreted  as  a  variable. 

The  language  uses  arithmatic  operators  and  a  set  of  functions 
to  perform  its  tasks.  MIMIC  uses  the  same  operators  as  FORTRAN, 

+,  *,  and  /  for  addition,  subtraction,  multiplication,  and 

division.  Note,  however,  that  the  exponential  operator  "**" 
does  not  exist  in  MIMIC.  Functions  are  specified  by  three  letter 
mnemonic  codes  (e.g.,  "SQR"  for  SQUARE  ROOT,  "INT"  for  INTEGRAL). 
Combining  these  in  proper  forms  allows  the  user  to  program  and 
solve  any  problem. 

As  a  time  simulation  language,  MIMIC  allows  only  one  independ¬ 
ent  variable,  time,  which  increments  automatically.  All  variables 
change  with  time  as  expected  from  an  analog  simulator,  and  this 
makes  integration  fast  and  simple. 

MIMIC  accepts  data  only  at  the  beginning  of  the  program. 

The  FORTRAN  read  statement  has  no  counterpart  in  MIMIC.  Thus, 

MIMIC  cannot  be  run  interactively,  all  jobs  run  in  "BATCH"  mode. 
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The  user  can  initialize  as  many  variables  as  desired,  but  all  should 
vary  with  time  in  course  of  execution.  MIMIC  accepts  constant 
terms  which  have  the  same  value  at  the  start  of  each  run,  parame¬ 
ter  terms  which  change  value  for  each  run,  and  arbitrary  functions 
of  one  or  two  variables  which  can  be  constant  or  parameter.  Data 
has  a  fixed  format. 

A  MIMIC  program  gives  output,  in  fixed  format,  at  specified 
intervals  of  the  independent  variable,  time.  Any  variable  in  the 
program  is  available  for  output,  as  well  as  any  arithmatic  expres¬ 
sion  . 

MIMIC  also  has  the  capability  of  generating  up  to  10  plots  per 
run  of  as  many  as  five  variables  against  any  other  variable, 
including  time.  The  MIMIC  processor  selects  the  scales  and  labels 
the  plots. 

MIMIC'S  most  useful  feature  is  its  centralized  integration 
function.  It  eliminates  the  burden  of  building,  coding,  and 
debugging  an  integration  scheme,  which  may  or  may  not  give 
accurate  results.  MIMIC  will  take  the  time  integral,  or  deriva¬ 
tive,  of  any  variable  in  the  program.  Systems  of  differential 
equations  lend  themselves  well  to  this  programming  tool. 

It  takes  very  little  time  to  become  proficient  with  MIMIC. 
Reference  7  gives  a  complete  description  of  the  language.  It 
explores,  in-depth,  MIMIC'S  workings,  and  describes  each  of  the 
functions  available  to  the  programmer.  This  powerful  tool  provides 
a  complete  time  history  simulation  of  any  physical  phenomenon 
quickly  and  easily. 
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APPENDIX  B  -  EXAMPLE  PROBLEM 


This  appendix  contains  a  complete  sample  run  of  the 
program  ENGOUT.  The  program  itself  is  listed,  as  is  all 
the  input  data  and  the  output  generated  by  the  program. 

The  example  is  run  for  a  theoretical  aircraft  with  one 
pair  of  main  gear  bogies.  It  is  powered  by  four  engines. 

The  spindown  time  constant  of  a  failed  engine  is  one-half 
second.  The  run  begins  at  zero  forward  speed,  and  full 
thrust  is  entered  as  a  step  function  at  time  zero.  The 
aircraft  is  equipped  with  a  yaw  damper  which  enters  rudder 
as  a  function  of  yaw  rate.  The  damper  has  a  maximum  authority 
of  four  degrees.  The  runway  is  dry  and  hard,  as  reflected  by 
the  Runway  Condition  Reading  (RCR)  of  23.  The  simulated  pilot 
has  a  reaction  time  of  one  second,  afterwhich  he  commands  the 
full  rudder  deflection  available  for  his  instantaneous  speed. 

The  aircraft  accelerates  to  a  speed  of  140  knots  true  airspeed, 
at  which  a  failure  of  the  outer  right  engine  occurs.  The 
center  of  gravity,  for  the  319,700  pound  gross  weight,  is  5% 
of  the  mean  aerodynamic  chord  forward  of  the  reference ,  as 
reflected  by  the  parameter  CG.  There  is  no  crosswind. 

The  input  data  appears  as  it  actually  would  on  a  printout. 

It  is  not  inputted  in  this  fashion,  but  in  the  MIMIC  format 
as  described  in  Appendix  A. 

Prior  to  failure,  output  is  printed  every  second  as 
specified  by  the  parameter  DTI.  The  output  shown  describes 
the  aircraft  accelerating  with  essentially  no  deviation  from 
the  runway  centerline.  The  variable  TRIM  remains  equal  to 
one  until  an  engine  failure  occurs,  at  which  time  it  is  set 
equal  to  zero,  36.43  seconds  after  starting. 

After  failure,  output  is  printed  every  0.2  seconds,  as 
specified  by  the  input  parameter  DT2 .  This  allows  closer 
tracking  of  parameters  of  interest.  For  a  full  second  after 
failure,  the  pilot  does  not  enter  directional  control,  and 
the  failure  induced  yaw  motion  is  opposed  only  by  the  aircraft 
directional  aerodynamics  and  the  yaw  damper.  The  damper  enters 
its  full  authority  to  combat  the  yaw  rate  induced  by  the  engine 
failure.  After  one  second,  the  pilot  commands  the  full  actuator 
output  to  enter  full  rudder  control.  The  rudder  input,  as  well 
as  nose  wheel  steering  commanded  by  the  pilot  thru  the  rudder 
pedals,  causes  a  yaw  rate  which  allows  the  pilot  to  regain 
control  of  the  aircraft,  and  eventually  realigns  the  aircraft 
with  the  runway  centerline  6  seconds  after  the  failure  occurred. 
The  aircraft  has  deviated  13.7  feet  from  the  centerline.  The 
last  action  of  the  program  is  to  print  the  final  values  of  all 
the  variables  in  the  program. 
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This  example  illustrates  the  use  of  ENGOUT.  It  is  of 
interest  that  the  failure-control  process  takes  only  6  seconds. 
This  rapidly  occurring  process  is  easily  affected  by  a  small 
change  in  almost  any  of  the  input  parameters. 
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***  MIMIC  SOURCE  RROCRAK  *• 
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